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ABSTRACT 

The origin and structure of magnetic fields in Gamma-Ray Burst (GRB) fireball plasmas are two of 
the most important open questions in all GRB models. Recent measurements of 7— ray polarization 
were claimed to suggest the presence of a uniform field originating in the compact object driving 
the outflow. This interpretation is, however, controversial, since a high degree of linear polarization 
is also possible in the presence of a random magnetic field, arguably originating in electromagnetic 
instabilities that develop at the collisionless shock. We show that the structure and strength of the 
magnetic field may be constrained by radio and IR observations of the early afterglow, where plasma 
effects on the polarization of the propagating radiation are significant. We calculate these propagation 
effects for cold and relativistic plasmas, and find that in the presence of a uniform equipartition field the 
degree of linear polarization is suppressed and circular polarization prevails at low frequencies, below 
~ 1 — 3 GHz (~ few x 10 14 Hz) in the forward (reverse) shock, at the onset of fireball deceleration. At 
higher frequencies linear polarization dominates. At the frequency of the transition between circular 
and linear polarization, the net level of polarization is minimal, ~ 10-20%. These features are nearly 
independent of the density of the environment into which the fireball expands. When the uniform 
field is much weaker than the equipartition value, the transition frequency is smaller by an order 
of magnitude. Depending on the geometry of the emitting region, oscillations of the polarization 
position angle may be observed from the optical reverse shock emission, provided that the strength 
of the magnetic field is close to equipartition. The dependence of these results on viewing geometry, 
outflow collimation and magnetic field orientation is discussed. When the field is entangled over 
length scales much smaller than the extent of the emitting plasma, the aforementioned effects should 
not be observed, and a linear polarization at the few % level is expected. Polarimetric observations 
during the early afterglow, and particularly of the reverse shock emission, may therefore place strong 
constraints on the structure and strength of the magnetic field within the fireball plasma. 
Subject headings: gamma rays: bursts — magnetic fields — plasmas — polarization — shock waves 



1. INTRODUCTION 

The past decade has been marked by major advances 
in our understanding of Gamma- Ray Bursts (GRBs). 
These bright explosions at cosmological distances, are 
now commonly believed to be generated by the dissipa- 
tion of a relativistic outflow associated with the gravi- 
tational collapse of a stellar mass object to a compact 
remnant (for reviews, see Piran 2000; Meszaros 2002; 
Waxman 2003a). Observations suggest that the out- 
flow is probably collimated to within a narrow jet, with 
smaller opening angles corresponding to apparently more 
energetic bursts (Frail et al. 2001). Current models re- 
late the production of 7— rays to emission of synchrotron 
or inverse-Compton radiation by relativistic electrons. 
Subsequently, the interaction of the ejecta with the cir- 
cumburst gas drives a forward shock into the ambient 
medium, and a reverse shock into the ejecta. This stage 
during which the forward and reverse shocks coexist is 
usually dubbed the "transition phase" ; once the reverse 
shock fully crosses the ejecta, the dynamics of the ex- 
panding fireball is adequately described by a selfsimi- 
lar solution (Blandford & McKee 1976). During the 
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self-similar phase which follows, a forward shock prop- 
agates in the ambient gas, gradually decelerating, un- 
til the inverse of the shock Lorentz factor drops below 
the jet opening angle (indicating the "jet break time"), 
at which stage the jet quickly slows down and spreads 
sideways. Subsequently the expansion continues sub- 
relativistically. Electrons accelerated at these shocks 
emit synchrotron radiation, which appears as an opti- 
cal flash from the reverse shock (Meszaros, Rees & Pap- 
athanassiou 1994; Akerlof et al. 1999; Zhang, Kobayashi 
& Meszaros 2003, and reference therein), lasting typically 
for the duration of the GRB itself, as well as a long-lived 
transient afterglow from the forward shock, characterized 
by a decreasing peak frequency, which lasts for weeks or 
even months (Frail, Waxman & Kulkarni 2000). 

Despite the successes of the models in predicting and 
explaining light curves and fluxes at various epochs and 
various wavebands, key questions, concerning predom- 
inantly the nature of the progenitor, the launching of 
the relativistic outflow, and the physics of the relativistic 
shock waves, are still unanswered. In addition, the physi- 
cal mechanisms responsible for the production of 7— rays 
are still subject to controversy. Leading scenarios for the 
production of the prompt 7— ray emission can be grouped 
into three main categories. The first consists of models 
in which the relativistic outflow is dominated by the ki- 
netic energy of the ejecta, and the 7— rays are the conse- 
quence of synchrotron emission by relativistic electrons, 
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which are Fermi-accelerated in optically-thin collisionless 
shocks created by inhomogeneities in the ejecta ("inter- 
nal shocks"). The assumption that the 7— rays are syn- 
chrotron photons requires strong magnetic fields, which 
probably result from the amplification of a small scale 
field by plasma instabilities in collisionless shocks (Gruzi- 
nov & Waxman 1999; Medvedev & Loeb 1999). In this 
scenario, the coherence scale of the magnetic field may 
grow to the size of the causally connected regions in the 
outflow, a size which is still smaller than the scale of the 
emitting slab. Afterglow observations strongly favor the 
formation of strong magnetic fields at relativistic shock 
fronts, and may therefore indicate that such a mechanism 
is also governing the emission from the internal shocks. 
In a second type of models the outflow is Poynting-flux 
dominated, and current-driven instabilities at large radii 
lead to acceleration of pairs that emit 7— rays by syn- 
chrotron radiation (Lyutikov, Pariev & Blandford 2003). 
The third type of models associate the prompt emission 
with the relativistic boosting of an ambient photon field, 
possibly related to an earlier SN explosion (Lazzati et 
al. 2000), by electrons outflowing with a high Lorentz 
factor. 

The recent detection of a strong linear polarization sig- 
nal in the prompt 7— ray emission of GRB 021206 by 
Coburn & Boggs (2003) triggered much interest in GRB 
polarization features [but see Rutledge & Fox (2003), 
who critically reanalyze the lightcurve of GRB 021206, 
and find no polarization signal]. This measurement - 
if real - provided a strong indication that synchrotron 
emission is indeed the underlying radiation mechanism 
of the GRB phenomenon (Waxman 2003b). Its interpre- 
tation in terms of magnetic field structure is, however, 
ambiguous. The high polarization in 7-rays was claimed 
to suggest the presence of a uniform magnetic field in the 
emitting plasma (Lyutikov et al. 2003), thus supporting 
the paradigm of a Poynting-flux dominated outflow ad- 
verting a "primordial" field from the "inner engine" . It 
turns out, however, that this line of reasoning is not free 
of caveats. First, it was recently demonstrated (Waxman 
2003b; Granot 2003; Nakar, Piran & Waxman 2003) that 
an outflow which is collimated to within a narrow jet and 
is observed off axis, may exhibit a high degree of linear 
polarization, even in the presence of a random magnetic 
field. Second, even the presence of a uniform field is 
still not an indication that the outflow is Poynting-flux 
dominated, since the amplitude of a toroidal magnetic 
field decreases as i?" 1 and tends to dominate at large 
radii - even if the outflow is dominated by kinetic en- 
ergy. The uniform field hypothesis is controversial also 
since it appears difficult to reconcile the existence of a 
large scale magnetic field with the presence of relativis- 
tic charge particles, which are accelerated in situ, since 
current understanding of particle acceleration requires 
the dissipation of these large-scale fields. 

In this paper we examine the fireball scenario, in which 
the outflow is dominated by the plasma kinetic energy, 
and demonstrate that the strength and structure of the 
magnetic field may be constrained by radio and IR ob- 
servations of the early afterglow, where plasma effects 
on the propagation of synchrotron radiation are signifi- 
cant. Specifically, we show that the presence of a uniform 
field permeating the plasma has distinctive fingerprints 
on the polarization properties of emitted radiation, which 



discriminate it from the case of a field entangled over 
small length scales. Of particular interest is the effect on 
the polarization of the reverse shock emission, since that 
plasma is just the ejecta producing the prompt 7— ray 
emission. Hence constraining the structure and strength 
of the magnetic field in the reverse shock is complemen- 
tary to measurements of the 7— ray polarization, and has 
direct implications for the plasma conditions in the ejecta 
and consequently for processes taking place in the vicin- 
ity of the compact source. This powerful probe may help 
to put constraints on models of field generation, thus pro- 
moting our understanding of both the shock physics and 
the nature of the "inner engine" . 

This article is organized as follows. In §[21 wc provide 
a brief review of polarized light transfer in a magne- 
toactive plasma, along with new results for a relativistic 
plasma. Next (§|3J) follows a qualitative discussion of fire- 
ball geometry and observation geometry, in the context 
of plasma propagation effects. The propagation effects 
are applied to GRB plasmas in §0J In S I4.1I we obtain 
plasma conditions during the early afterglow, consider- 
ing both cases of an expansion into a uniform-density 
interstellar medium (ISM) in § 14.1.11 and an expansion 
into a circumburst wind in § 14.1.21 Detailed analysis of 
propagation effects at the corresponding shocks follows 
in § 14.21 The implications of the strength and structure 
of the magnetic field are discussed in § 14.41 Finally, we 
consider the implications of our results in §[S] 

2. PLASMA EFFECTS ON THE POLARIZATION OF 
PROPAGATING SYNCHROTRON RADIATION 

As is well known, a plasma permeated by a magnetic 
field is birefringent, and a phase shift accumulates be- 
tween the normal modes as radiation propagates through 
the plasma. The observational outcome of this phe- 
nomenon is the Faraday rotation effect. Below we briefly 
review the formalism of transfer of polarized synchrotron 
radiation in a magnetoactive plasma f§ l2.1|) . and bring 
new results for a plasma dominated by a relativistic elec- 
tron population f§ l2.2fl . The appendix gives a more com- 
prehensive derivation of these results. 

2.1. Formalism 

The effect of the magnetoactive plasma on the transfer 
of polarized radiation may be treated most straightfor- 
wardly via the analysis of the wave equation, incorporat- 
ing the effect of the medium by a dielectric tensor. When 
the radiation can be regarded as transverse, the tensor 
K ij {hi = 1,2) 4 representing the plasma effect in the 
dielectric tensor, is separated into a sum of hermitian 
and anti-hermitian tensors as follows : 

A)- m 

The anti-hermitian component of the tensor is responsi- 
ble for absorption. The hermitian part is responsible for 
propagation effects; specifically, / = Afc(cj) = wAn(w)/c 
is the difference between the wave numbers of the nor- 
mal modes, giving rise to the Faraday rotation effect. In 
the case of synchrotron self-absorption, the parameter v 
is typically smaller than k by a factor 7 m , where "f m is 

4 We adopt a right-handed system of coordinates with the wave- 
vector k directed along axis 3. 
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a typical Lorentz factor of the relativistic electrons, and 
will therefore be neglected hereafter. 

The wave equation for the Fourier components of the 
electric field can be manipulated, with the help of equa- 
tion Q, into a transfer equation for the Stokes parame- 
ters : 

dJ/ds = Ei — nl — qQ 
dQ/ds = e Q -qI- K Q-fU 
dU/ds = fQ-nU-hV 
„ dV/ds = e v + hU - kV , 

where s is a length parameter along the path of the ray, 
and ej,eq and ev are the emission coefficients corre- 
sponding to /, Q and V. 

The physical significance of the propagation parame- 
ters / and h is most easily demonstrated in a transpar- 
ent medium. As an example, consider the case f h, 
with negligible k and q. Equation JQl then implies that 
the normal modes are (left- and right- ) circularly polar- 
ized. Integrating equation J5J we obtain oscillations in 
the Stokes parameters Q and U, leading to oscillations 
and damping in the degree of emergent linear polariza- 
tion : 
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A large phase shift accumulated between the normal 
modes, A(j) = /s > 1, results in damping of the de- 
gree of linear polarization by a factor of A(f>. This is 
accompanied by oscillations of the polarization position 
angle x — I tan _1 ([//(5), providing the Faraday rotation 
effect. Analogous, simple to interpret phenomena char- 
acterize also the case h 3> / (linearly polarized normal 
modes). 

2.2. New results for relativistic plasma 

The coefficients / and h are easily derived when the 
plasma is cold : 
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2clu(lo 2 — uj 2 b ) 

where oj 2 = 4-7m e e 2 /m e is the non-relativistic elec- 
tron plasma frequency, and ujb — eB/m e c is the non- 
relativistic electron Larmor frequency -d is the angle 
between k and B, measured in the frame comoving with 
the plasma. 

In the context of GRB plasmas, however, one also en- 
counters a situation in which the plasma is dominated 
by a relativistic electron population. In our subsequent 
derivation of analytic expressions for the propagation co- 
efficients in a relativistic plasma we assume that : (i) the 
electrons are highly relativistic, characterized by a spec- 
trum with a power-law tail extending to high energies; 
(ii) their distribution function is isotropic; and (Hi) the 
deviation of the dielectric tensor from its vacuum value 
is small. The first assumption is implied by afterglow 
observations, whereas the second is adopted for simplic- 
ity. The third assumption will be justified later for the 
specific parameters of GRB plasmas. In addition, we 
restrict the discussion to a frequency range u> ^> ll>b, 



so that modifications to synchrotron emission and self- 
absorption coefficients due to the discrete nature of low 
harmonics of vb could be safely ignored. 

Our main results are summarized in the following equa- 
tions (for an outline of the derivation, see Appendix 
E3J : 
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Here v = lo/2it, j m is a characteristic electron Lorentz 
factor (typically that of the thermal electrons), v m — 
(3eB) I (Airmec)^ is the characteristic synchrotron fre- 
quency corresponding to j m , and p is the spectral in- 
dex of the power-law tail of the electron distribution 
function. Equation JHJ) generalizes an earlier result by 
Sazonov (1969), who considered only the frequency range 
v > v m . 

After the completion of this work, we learned of an in- 
dependent derivation of similar expressions for the prop- 
agation coefficients in the relativistic regime by Mat- 
sumiya and Ioka (2003). Our expressions [eqs. © and 
agree nicely with theirs; nevertheless the resulting 
polarization differs substantially, as we further discuss in 
§0 

Using these results in the context of Faraday rotation, 
we notice that /, and therefore also A</>, are proportional 
to v~ 2 . The frequency interval A^, corresponding to a 
single period of the position angle oscillation, translates 
to an increase of 2n in A</>, from which we get (for a cold 
plasma) : 

a, ....^v 

e 3 n e BW ' V ' 

where W is the width of the slab through which the ra- 
diation propagates. The generalization for a relativistic 
plasma is trivial. 

3. CONSEQUENCES OF FIREBALL GEOMETRY AND 
VIEWING GEOMETRY 

A strong polarization signal clearly requires a depar- 
ture from spherical symmetry, either via a non-spherical 
outflow configuration and a special location of the ob- 
server with respect to it, or through a preferred direction 
of the magnetic field in the plasma (or, of course, a com- 
bination of the two). Geometrical considerations related 
to the fireball collimation, viewing angle and field's direc- 
tion are essential for the calculations we describe below, 
and so before presenting detailed results of these calcu- 
lations, we make a short digression, and discuss qualita- 
tively some important consequences of the geometry. In 
particular, we shall show that the angle on which the 
propagation coefficients / and h depend [c.f. equations 
01 - ©] is not likely to assume special values such as 
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and 7r/2 in a typical GRB jet, and may approach only 
for a narrow jet observed off-axis, with a magnetic field 
oriented normal to the shock plane; These two regimes 
have markedly different observational consequences. We 
focus in this section only on the case of a large scale, 
uniform magnetic field. The case of a fluctuating field is 
treated in H4.4I 

Fig. ^ shows schematically the essential ingredients 
of the fireball and viewing geometry. The figure also 
illustrates (with shaded areas) electron cooling due to 
synchrotron losses away from the (forward and reverse) 
shock fronts. This effect, which has important conse- 
quences for the resulting polarization, is discussed in 
greater detail in § 14.21 
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Fig. 1. — Schematic illustration of the fireball and viewing ge- 
ometry. 9j and 9 a are, respectively, the jet opening angle and the 
angle between the line of sight to the observer and the jet axis. Ti 
is the Lorentz factor of the ejecta prior to deceleration, whereas T 
is the Lorentz factor of the plasma shocked by the reverse and for- 
ward shocks. The initials RS, CD and FS denote the reverse shock, 
the contact discontinuity and the forward shock, respectively, while 
CBM denotes the circumburst material. The dark shading marks 
the cooling of the electrons, with darker shade representing colder 
electrons. 



Mounting evidence strongly suggests that the relativis- 
tic outflow in GRBs is collimated into a jet, the opening 
angle of which is usually inferred based on estimates of 
the jet break time from optical observations (Frail et al. 
2001). These observations are carried on a day timescale 
after the GRB trigger, when the jet has already decel- 
erated to a low Lorentz factor, T ~ 10, whence one de- 
duces 9j ~ 0.1. We, on the other hand, are interested 
in the emission from the early afterglow, at which stage 
the fireball expands with a much higher Lorentz factor, 
r ~ 10 25 [see eqs. (|SJ) and <|14[) ]. R is unclear whether 
during this early stage the jet opening angle is similar 
to the one inferred from optical observations. However, 
Freedman & Waxman (2001) have demonstrated that the 
isotropic-equivalent fireball energy may be robustly esti- 
mated from the X-ray afterglow flux on a day timescale 
after the GRB, and furthermore that this estimate ex- 
hibits a strong correlation with the isotropic-equivalent 
energy inferred from the GRB 7-ray fluence [see also 
Berger, Kulkarni & Frail (2003) for a recent analysis 
compiling more X-ray data]. This correlation, and in 
particular the fact that the latter estimate is never much 
larger than the former, strongly argue in favor of a sim- 
ilar jet opening angle during the GRB and early after- 



glow phases. Since the X-ray afterglow is obtained on 
a day timescale, when the jet has already decelerated 
to T ~ 10, this arguments suggests an opening angle 
9j ~ 0.1 for a typical GRB jet. An observer of the early 
afterglow is therefore not likely to see the edge of a typi- 
cal jet, and is consequently not sensitive to deviations of 
the fireball from spherical symmetry (other than a large- 
scale magnetic field in a particular direction). 

Deviations from spherical symmetry may be mani- 
fested for an observer located on a line of sight which 
makes an angle ~ T -1 from the jet edge, and are there- 
fore probable only if the jet is extremely narrow, with 
0j ~ r _1 . During the GRB and early afterglow phases 
the fireball expands with T ~ 10 25 , and so the argu- 
ments presented above suggest that such narrow jets are 
not generic; However, in the framework of the "stan- 
dard energy reservoir" hypothesis promoted by Frail et 
al. (2001), brighter GRBs are naturally associated with 
narrower jets. Adopting this energy-angle relation, GRB 
021206, which was among the strongest bursts ever, had 
an opening angle < 1/40 rad, or even smaller (Nakar et 
al. 2003). This observation is key in associating the high 
degree of linear polarization measured in GRB 021206 
with an observational alignment effect (Waxman 2003b, 
Nakar et al. 2003). 

Upon transforming to the frame comoving with the 
fireball, the distinct geometries of a typical jet and a 
narrow jet correspond to the radiation k vector making 
either a large T _1 ) or a small (~ T _1 ) angle relative 
to the plane of the shock, for a typical (9j ^> T _1 ) and 
narrow (8j ~ T _1 ) jets, respectively. Assume now that 
the magnetic field is dominated by a component either 
parallel or normal to the plane of the shock. The angle 
$ between k and B may assume special values, i.e. close 
to or to 7r/2 (see an illustration of these two geome- 
tries in Fig. |2J). Recalling the angular dependencies of 
the propagation coefficients, / oc cos 1) and h cx sin 2 1? 
[approximately, see equation ©] , this may result in very 
different propagation effects, and hence in qualitatively 
different polarization properties of the emitted radiation, 
as we explain next. 





Fig. 2. — Geometrical setup in the frame comoving with the 
plasma. Shown are k — the direction of the radiation wave vector 
(i.e. the line of sight), — a component of the magnetic field 
normal to the shock plane, and By — the field component lying 
in the shock plane. The angle $ between the line of sight and the 
field (either component) is also shown. Left : A typical viewing 
geometry for the case 9j S> r — 1 ; The wave vector k lies outside 
the plane of the shock. Right : The viewing geometry for the case 
9j ~ r~ ; Now k lies roughly in the shock plane. 

In § 14.21 we show that the frequency ranges where ob- 
servational fingerprints of plasma propagation effects are 
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expected are v ~ 1 GHz for the forward shock emis- 
sion, and v ~ 10 13 Hz for the reverse shock emis- 
sion. Let us assume that the plasma is permeated by 
a large scale magnetic field, and denote by Rf/h — 
(f/h) sin 2 fl/cosd the (geometry-independent) ratio of 
the propagation coefficients in the relevant frequency 
ranges. The plasma parameters we derive in § 14.21 im- 
ply that Rf/h ~ lOfV/lGHz) -1 / 3 in the forward shock, 
and Rf/ h r~j 10 3 (V/10 13 Hz) in the reverse shock 5 . Con- 
sider now a typical jet with Oj The high Rf/h 
imply that whatever the direction of the magnetic field 
be, in order for h to dominate over / (i.e. in order for 
■& to be sufficiently close to 7r/2), k must lie within a 
very narrow solid angle. Consequently, a configuration 
where / 3> h is much more probable than one in which 
h > / [by a factor of ~ R f/h - 1 ~ 9 (10 3 ) in the 
forward (reverse) shock]. Next, consider a narrow jet 
observed slightly off-axis, with 9j ~ Recall that 
in this regime, the vector k is lying in the shock plane 
(see Fig. [21 right). If the uniform magnetic field also 
lies predominantly in the shock plane, we find again that 
/ ^> h is much more probable than h ^> / [by factors of 
~ 7 and ~ 10 25 for the forward and reverse shocks, re- 
spectively]. If, however, B is normal to the shock plane, 
it is necessarily normal also to k, and therefore h ^S> f. 
The normal modes of propagating radiation are then lin- 
early polarized, and the resulting emission has markedly 
different polarization properties. 

In the next section we calculate propagation effects in 
both geometric regimes. The case / >■ h is studies in 
more detail, since this is the characteristic case for a typ- 
ical GRB jet. The polarization pattern for this regime 
is shown in Fig. In Fig. we show the polariza- 
tion properties of the h 3> / regime, which, as we argued 
above, is probable only for a narrow jet observed off-axis, 
with a uniform magnetic field normal to the shock plane. 
The qualitative results of this discussion, complemented 
by some quantitative consequences of the actual calcula- 
tions, are summarized in Tabled 

4. APPLICATION TO GRB AFTERGLOWS 

The application of the results derived in §|5] to plas- 
mas of GRB afterglows requires knowledge of the plasma 
number density, strength and structure of the magnetic 
field, and the electron distribution function produced by 
the shock. These parameters, in turn, depend on the 
total energy E carried by the relativistic plasma, the 
number density n of the ambient gas, the duration T 
of the burst, the Lorentz factor Ti reached by the out- 
flow prior to the deceleration phase, the fractions e e and 
€b of the energy density carried by relativistic electrons 
and magnetic field, respectively, and the spectral index 
p = — dlnn e /dln7 e of the assumed power- law distribu- 
tion of the shock-accelerated electrons. When consider- 
ing a fireball expanding into a wind-stratified medium, 

5 The reason for the much larger value of Rf/h m the reverse 
shock is essentially the fact that reverse shock radiation travels 
through a thick slab of cold plasma, whence the values of / and 
h follow from equation rather than from eqs. J2J and JfjJ, as 
they do in the forward shock. The values we give here for Rf/h 
are obtained from plasma parameters of a fireball expanding into a 
uniform-density ISM. Nevertheless, the corresponding values for a 
fireball expanding into a wind are similar, and hence this scenario 
will not be treated separately in this section. 



n is replaced by M /4irm p R 2 v w , involving the ratio be- 
tween the progenitor's mass loss rate, M, and the wind 
velocity, v w (where m p is the mean atomic mass for the 
gas). We shall therefore briefly address the observational 
constraints on these parameters, on which our choices of 
parameter values is based. 

We begin with the parameters characterizing the ex- 
plosion. First, we focus on the subclass of long-duration 
(T > 2 s) GRBs, since these are the only bursts localized 
by the BeppoSAX and HETE-II missions with identified 
afterglows. We therefore express our results in units of 
Tj = T/IO s. 

Second, as we are interested in the emission prior to 
the jet break time, the energy scale determining the dy- 
namics is not the actual explosion energy, but rather 
its isotropic equivalent, i.e. the total energy derived as- 
suming isotropic emission. Examination of the list of all 
bursts with measured redshifts indicates E 1 ,- lso ~ 3 x 10 53 
ergs as the mean isotropic equivalent energy inferred 
from the 7— ray fluence (Bloom, Frail & Kulkarni 2003). 
Since this value must be smaller than the total kinetic 
energy carried by the plasma, it is reasonable to express 
E- lso in units of 10 54 ergs. Hereafter we drop the subscript 
"iso" for brevity. 

Considering now the ambient gas density, when the ob- 
served light curve can be reconciled with expansion into 
a uniform density ISM, Bloom et al. (2003) find n to 
typically lie in the range 10°' 5±1 cm -3 . We shall there- 
fore adopt a parametrization with no = n/(l cm -3 ). On 
the other hand, it is now commonly believed that the 
long duration GRBs are probably associated with core- 
collapse supernovae. This association was originally mo- 
tivated by evidence of intense star formation in iden- 
tified GRB host galaxies (Bloom, Kulkarni & Djorgov- 
ski 2002), and by additional evidence for optical super- 
novae emission in several GRB afterglows (Bloom 2003); 
finally, the SN-GRB relation was put on solid grounds 
with the recent detection of a supernova emission spec- 
trum (SN2003dh) associated with GRB030329 (Stanek 
et al. 2003; Hjorth et al. 2003). Considering a likely 
phase of mass loss during the late stages of the mas- 
sive star's life, it is natural to examine a scenario in 
which the fireball expands into a wind medium. Below 
we consider the simplest wind model, where prior to the 
explosion the progenitor star ejects mass at a constant 
rate M = 10 -5 M-sM^yr -1 and at a constant speed 
v w = 10 3 w 3 kms _1 (Chevalier & Li 1999). These pa- 
rameters agree with observations of massive stars which 
are believed to be the progenitors of SNe Ib/c associated 
with GRBs (Willis 1991). Note, that since n ~ 1 cm" 3 
for a uniform-density circumburst gas is inferred from 
observations carried typically ~ 10 hours after the burst, 
this is consistent also with a wind environment, which 
should have comparable density at the radius reached by 
the fireball on a ~ 10 hours time scale (Livio & Wax- 
man 2000). Finally, the Lorentz factor reached by the 
outflow prior to its deceleration by the circumburst gas 
is typically a few hundreds. We choose to parameterize 
our results in units of 1^2.5 = I^/IO 2 ' 5 . 

We now turn to consider the parameters p, e e and 6b- 
There is strong indication that the p and e e are uniform 
among different bursts, with e e close to an equipartition 
value (Frail et al. 2001; Freedman & Waxman 2001; 
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TABLE 1 

Effect of fireball, observing and field geometry on the polarization 





Uniform field 




Random field 


_L a 




II b 


_L | 



8j » r- 1 / > h f>h |/| > \h\ 

Faraday depolarization Faraday depolarization no propagation effects 

FS c : n c > U L at v < 1 GHz U c > Tl L at u < 1 GHz no polarization n < 10% 

RS d : U c > U L at v < few X f 14 Hz Yl c > U L at i/ < few X 10 14 Hz no polarization n < 10% 



/>»/ f>h \h\>\f\ \f\>\h\ 

no Faraday depolarization Faraday depolarization no propagation effects 

FS : n c - n L - 20% - 50% at v < 1 GHz n c > Tl L at v < 1 GHz n c < n L - n < few x 10% 

RS : n c < n L < 100% at v < few x 10 13 Hz U c > n L at f < few x 10 14 Hz n c < n L ~ n < few x 10% 



a Ficld normal to plane of shock 
a Ficld parallel to plane of shock 
c Forward shock emission 
^Reverse shock emission 



Berger, Kulkarni & Frail 2003), and p = 2.2 ± 0.1 (Wax- 
man 1997a; Galama et al. 1998; Stanek et al. 1999; Frail, 
Waxman & Kulkarni 2000; Freedman & Waxman 2001). 
We therefore adopt e e — 0.1e e ,-i and p — 2.2 in the nu- 
merical examples we consider below. The value of sb is 
less constrained by observations, and its estimated value 
ranges from e B ~ 10 _1 (e.g., Waxman 1997a; Wijers 
& Galama 1999) to e B - 10~ 6 (e.g., Wijers & Galama 
1999; Chevalier & Li 1999; Galama et al. 1999; Wax- 
man & Loeb 1999). Since in cases where es can be 
reliably constrained by multi-waveband spectra, values 
close to equipartition are inferred (Frail et al. 2000), 
we first explore consequences of propagation effects on 
emission from plasmas with eg = O-le^.-i- This is also 
the parametrization used in S I4.1I We then study the 
implications of a weak magnetic field in § 14.41 adopting 
€b = 10~ 4 as an illustrative example. 

Next we explore how the propagation effects described 
in the previous section manifest themselves in the con- 
text of a GRB afterglow. In S I4.1I we derive the relevant 
plasma conditions at the onset of fireball deceleration, 
considering both expansion into a uniform density ISM 
fS I4.1.1[) and into a wind medium f 5 14.1. 2f> . In 5 14.21 we 
examine the observable consequences of our derivations, 
investigating first the case of a uniform, close to equipar- 
tition magnetic field. Since a typical GRB jet is charac- 
terized by 9j ^> r _1 , implying typically / ^> h in the 
frequency ranges of interest (see JjEJ, this section gives 
detailed analysis of this regime. 5 l4.3l studies the comple- 
mentary case, where h 3> /. The effect of a weak field, 
as well as the effect of an entangled field geometry, are 
treated in S 14.41 Throughout the rest of this paper, we 
use primed quantities to denote values measured in the 
frame that is comoving with the plasma. 

4.1. Plasma parameters during the early afterglow 

4.1.1. Expansion into a uniform- density ISM 

The relativistic blast wave driven into the ambient gas 
by the energetic explosion approaches a self-similar be- 
havior (Blandford & McKee 1976) once a reverse shock 
crosses the ejecta and heats it (Meszaros & Rees 1997). 



During this transition episode both the shocked ambient 
gas and the heated ejecta propagate with a Lorentz fac- 
tor which is close to that given by the Blandford-McKee 
self-similar solution (see, e.g. Waxman & Draine 1999) : 

r W ~ ( 17 E - K ~ 328 EH S n- 1/8 T7 3/s , 

\ 10247mm p c 5 T 3 / 54 1 ' 

(8) 

where E54 — E/10 5 erg is the isotropic-equivalent of the 
fireball energy, T\ — T/lOs is the burst duration, and 
no = n/lcm -3 is the ambient number density. The su- 
perscripts R, F denote the reverse and forward shocked 
plasmas, respectively. The corresponding electron num- 
ber densities are accordingly given by the shock jump 
conditions : 

n'W =4Tn ~ 1300 E^n^T^ 8 cm" 3 

(rVr<K^ 

~ 3.99 x 10 5 EH 8 nl /8 Ti g/S rri 5 cm - 3 . (9) 

The comoving widths of shocked forward and reverse 
plasma shells can be shown to be comparable, and are 
estimated as : 

W '(F,R) ^ R/AV „ 9 8 x 1Q 13 £l/8 n -l/8 r 5/8 cm ^ (1Q) 

where R ~ 1.3 x 10 17 cm is the radius where the reverse 
shock is propagating through the ejecta. From the shock 
jump conditions we infer the strength of the magnetic 
field : 

B'( F >V = [32ire B nm p c 2 ] 1/2 r 

~ 40 £ 5 T^o /8T r 3/8 4 /2 -i G • ( n ) 

The distribution function of relativistic electrons injected 
by the shocks is assumed to increase as 7g up to a charac- 
teristic Lorentz factor 7 m (Gruzinov & Waxman 1999), 
and decrease as a power law r )~ p (with p = 2.2) at higher 
values of the Lorentz factor. Then the jump conditions 
imply : 

(F) = 4(p-2) mp 
lm 3(p-l) e m e 
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-1.32 x 10 4 £^^ 1/8 T 1 ~ 3/8 e e ,-x 

'(-R) = 4 (P~ 2 ) m P r ' 
7m 3(p-l) ee m e r 

~ 43 ^ 1 / 8 nJ /8 T? /8 r i , 2 . 6 e e ,_ 1 . (12) 

The observed characteristic synchrotron frequencies 
corresponding to the Lorentz factor of the thermal elec- 
trons in these two shock scenarios are : 

„W> ~ 1.04 x 10 14 r^^e 2 Hz; (13) 

In the parlance of the high-energy and IR communities, 

(-F) 

respectively, these read an energy El h ' =39.7 KeV for 

( H) 

the forward shock, and a wavelength Am = 2.9/im for 
the reverse shock. 

4.1.2. Expansion into a wind 

We now examine a fireball expanding into a wind. 
The simple wind model considered here results in a non- 
homogeneous ambient gas with a profile n oc Rr 2 . The 
reverse shock crosses the ejecta at a typical radius of 
2.2 x 10 16 cm, at which time the plasmas shocked by 
the forward and reverse shocks propagate together with 
a Lorentz factor 

~ 135 ^i 4 (M_ 5 / U3 )- 1/4 T- 1/4 . (14) 

Following the reasoning employed above, we obtain the 
other relevant plasma parameters for the wind scenario : 
W >(F,R) ^4. 4 x 10 i3 J B 5 1 4 / 4 (Af_ 5 / U3 )-i/4 r 3/4 cm 

<< F >~3.43 x 10 5 £-/ /4 (Af_ 5 A, 3 ) 7/4 Tr 5/4 cm" 3 
n'W-1.78 x 10 7 E^iM^/vs^T^TT^ cm" 3 
B'^ -418 ^ 54 1/4 (M_ 5 /^ 3 ) 3/4 T- 3/4 e^ 2 _ 1 G 

7 ^^c5440i? 1 4 / 4 (M_ 5 M)- 1/4 rr 1/4 ee ,- 1 

</„[*> ~ 105 B^iM^/v^X^^e,-! ■ (15) 

The characteristic synchrotron frequencies in the wind 
scenario are : 

,Cf) * 7.0 x 10- E^T-^el-^T-i Hz - 

4?) ~ 2.60 x 10 15 S B - a 1 / a (M_ 6 M)T 1 - 1/2 l1 , ^. 6 

x^e^Hz. (16) 

4.2. Observational consequences for a uniform, close to 
equipartition magnetic field — the "typical jet" 
regime 

The formalism developed in S I2.1I was applied to the 
plasma conditions derived above. We assumed the mag- 
netic field and electron density are uniform across the 
slab. As an electron drifts away from the shock front 
with a comoving velocity c/3, it loses energy due to syn- 
chrotron emission. Consider, for simplicity, an isotropic, 
monoenergetic electron population, with a Lorentz fac- 
tor 7j' endowed by the shock. Assuming this population 
remains isotropic (in momentum space of the comoving 
plasma) as it loses energy to radiation, the mean Lorentz 
factor decreases as 



1 + 3<7 T £ /2 7^' /fam e c 2 ' 



where z' is the comoving separation from the shock front. 
Assuming now that the shock injects a spectrum of elec- 
tron energies, equation (|17l) implies that the shape of 
the energy distribution function is modified as the sepa- 
ration from the shock front increases. As a result, only a 
thin slab of width 3 x 10 11 cm behind the forward shock 
is dominated by a relativistic electron population, while 
the reverse shock front is separated from the observer by 
a cold slab 2 x 10 13 cm wide (for a fireball expanding 
into ISM; The corresponding widths for the wind sce- 
nario are 5 x 10 9 cm and 4 x 10 13 cm, respectively). As 
this strongly affects the values of emission, absorption 
and propagation coefficients (see Appendix lA")) . we in- 
tegrated the radiative transfer equation @ by dividing 
the emitting slab into a large number of sub-slabs, such 
that within each sub-slab the electron energy spectrum 
did not change significantly. The coefficients of emission, 
absorption and propagation could then be assumed con- 
stant within each sub-slab, allowing an analytic solution 
of the transfer equation. This division was refined until 
convergence was obtained. 

Equations Q, JSJ) and © imply a strong dependence 
of the propagation coefficients on the comoving angle 
between the directions of the wave-vector and the mag- 
netic field. Yet, as argued in §[21 the observing geometry 
for a typical GRB jet are likely to give rise to a situation 
where / ^> h. In such circumstances •& — 7r/2 is most 
unlikely, and the angular dependence of / and h may be 
ignored, for simplicity. This simplification was adopted 
in the calculation descried below, and gives a good or- 
der of magnitude estimate of the emitted polarization of 
a typical jet with probable viewing geometry. The case 
/i> /, occurring when the jet is narrow (0j ~ F -1 ) and 
the magnetic field is predominantly normal to the shock 
plane (see §01 is treated separately in § 14.31 

We verified that all elements of the tensor Kij at all 
relevant frequencies are much smaller than unity, thus 
self-consistently justifying the analysis of § l2.1l which as- 
sumed that the electric field is perpendicular to the wave 
vector. 

The resulting emitted linear, circular and total degrees 
of polarization at the onset of fireball deceleration are 
shown in Figs. (forward shock) & 0] (reverse shock). 
We show results for both ISM and wind scenarios. All 
four graphs exhibit broadly the same behavior. Due to 
the Faraday effect, the linear polarization of emitted ra- 
diation is strongly suppressed at low frequencies, and 
circular polarization dominates. At a frequency of — 1 
GHz (30 GHz) in the forward shock into ISM (wind), and 
— 3 x 10 13 Hz in the two reverse shocks, linear polariza- 
tion takes over. At the transition frequency the total 
degree of polarization has a global minimum (10-20%). 
Only at frequencies much higher than this transition fre- 
quency the canonical degree of polarization (dominated 
by its linear component) is reached. 

The emission from the reverse shock overtakes the for- 
ward shock emission at frequencies above v ~ 2 x 10 10 
Hz. We find that for the reverse shock, the transition 
from circular to linear polarization is accompanied by 
a strong oscillation of the polarization position angle as 
function of frequency. Due to synchrotron losses, most 
of the plasma heated by the reverse shock cools quickly 
on a dynamic time scale, and consequently Faraday rota- 
tion by the cold plasma strongly affects radiation emitted 



from regions closer to the reverse shock front. The re- 
sulting oscillation is a hallmark of reverse shock emission 
up to UV frequencies (10 15 — 10 16 Hz), and is character- 
ized by Av/u ~ lQ~ 1 v'l 5 [from equation Q], where v\$ is 
the observed frequency in units of 10 15 Hz; Observation 
of these oscillations may therefore be limited by instru- 
ment resolution in the optical and particularly the IR 
bands. Note, additionally, that during the transition to 
self-similar dynamics, emission is observed from an area 
on the sky with a radius comparable to R/T (Waxman 
1997b). Hence radiation at observed frequency v spans 
a range of comoving frequencies Az/ . Since radiation 
at different comoving frequencies is subject to different 
phase shifts between normal modes, spectro-polarimetric 
features characterized by a smaller comoving Az/ may be 
washed out altogether. Although an exact, quantitative 
statement regarding the implications to the observabil- 
ity of the Faraday oscillations is beyond the scope of this 
paper, we expect that the depolarization of the linear 
component, caused by either low spectral resolution or 
smearing of the signal by the source plasma, is expected 
to yield a signal dominated by a circular polarization 
He ~ 0.1(^/10 13 Hz) -1 / 3 , surviving up to UV frequen- 
cies. 

A superposition of the polarization patterns of the for- 
ward and reverse shocks is shown in Fig. [SJ for the ISM 
scenario. The polarization is dominated by the reverse 
shock emission at frequencies above ~2x 10 11 Hz. The 
analogous superposition of forward and reverse shock po- 
larization for the wind scenario is qualitatively similar. 

Remarkably, the main characteristics of the polariza- 
tion pattern, including the frequency of transition from 
circular to linear polarization, are largely independent of 
the circumburst density, and are dictated predominantly 
by the nature of the shock (i.e. whether the shock is for- 
ward or reverse). This is particularly intriguing in view 
of the three orders of magnitude factor separating the 
densities of ambient gas in the ISM and wind scenarios, 
at the onset of fireball deceleration. Considering the two 
forward shocks first, it is straightforward to show that 
in both cases v' a -C v' m , where v' a is the (comoving) syn- 
chrotron self-absorption frequency. At smaller frequen- 
cies, v' < z/, the effective width W' e ^{v') from which 

photons are emitted scales as W' e ^(v') oc [k,'{v')]~ x oc 
z/ 5 ^ 3 . By virtue of equations 10 and then, we have 
at low frequencies ILl(v') oc W' e ^v' 2 oc v' 1 ^ . On the 
other hand, the degree of circular polarization at frequen- 

2/3 

cies v' 3> v' B scales as v' 1 . Denoting the frequency of 
the transition from circular to linear polarization by v t 
we therefore find v' t oc (^^a) 1 ^ 2 - Transforming to the 
observer's frame, we find for the forward shock in the 
ISM, v t oc n 27 / 80 ~ n 1 / 3 . Substituting here the ambient 
density for the wind scenario at the transition-phase ra- 
dius, which is larger by ~ 3 orders of magnitudes, and 
correcting for the much narrower slab which effectively 
emits photons in this scenario (due to the more efficient 
synchrotron cooling), we obtain v t in the wind scenario 
larger by a factor of ~ 3, as can be readily verified in 
Figure |21 Considering now the two reverse shocks, the 
frequency at which linear polarization takes over circu- 
lar polarization is simply the synchrotron self-absorption 
frequency, since at v > v a the region closest to the shock 
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Fig. 3. — Propagation effects on synchrotron radiation propagat- 
ing in a magnetoactive plasma. The plasma conditions are charac- 
teristic of forward shocks for a fireball expanding into a uniform- 
density ISM (above) and into wind (below), during the onset of 
fireball deceleration. The magnetic field is assumed uniform and 
close to equipartition (eg = 10 -1 ). Both figures show the degrees 
of linear polarization II^ (solid), circular polarization Uq (dotted), 
and the total degree of polarization IT (dashed). 



front becomes observable. As this region is populated by 
electrons which are still highly relativistic, it contributes 
most of the photons above v a , while the thick slab of cold 
plasma separating this emitting region from the observer 
is characterized by a much smaller emissivity. Radiation 
emerging from the emitting region close to the reverse 
shock front is subject to a large Faraday rotation due to 
the presence of the cold plasma; however, since the emis- 
sivity of the cold plasma is small, this has only a minor 
effect on the degree of linear polarization of that radia- 
tion. The similar v t in the two reverse shock scenarios is 
thus a consequence of a similar self-absorption frequency 
in these two cases. 

We note that our calculations indicate the plasma 
shocked by the forward shock does not affect the polar- 
ization properties of the emission from the reverse shock, 
since it is optically thin in the frequency range where 
the reverse shock emission exhibits interesting polariza- 
tion features. Moreover, only the non-relativistic compo- 
nent of the plasma behind the forward shock introduces 
a Faraday rotation with a large enough amplitude above 
a frequency of 10 13 Hz. Yet this plasma contributes only 
negligibly to the emission and absorption of radiation, 
and hence does not change the linear, circular or total 
degrees of polarization, but merely modifies weakly the 
pattern of position angle oscillations. 

4.3. The "narrow jet" regime — observational 
consequences 

A narrow jet (6j ~ T _1 ) observed slightly off-axis, 
with a uniform magnetic field normal to the plane of 
the shock(s) is the only geometric configuration where 
the propagation coefficient h is very likely to dominate 
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Fig. 4. — Propagation effects on synchrotron radiation propagat- 
ing in a magnetoactive plasma. The plasma conditions are charac- 
teristic of reverse shocks for a fireball expanding into a uniform- 
density ISM (above) and into wind (below), during the onset of 
fireball deceleration. The magnetic field is assumed uniform and 
close to equipartition (eg = 10 -1 ). Both figures show the degrees 
of linear polarization II £, (solid), circular polarization Uq (dotted), 
and the total degree of polarization II (dashed). 




Fig. 5. — Polarization pattern of combined (superposed) re- 
verse and forward shock emission, for a fireball expanding into a 
uniform-density ISM. Shown are the degrees of linear polarization 
TIl (solid), circular polarization lie (dotted), and the total degree 
of polarization II (dashed). At frequencies u < 2 X 10 11 Hz the po- 
larization is determined by forward shock emission, whereas above 
that frequency it is dominated by the emission from the reverse 
shock. In the frequency range 3 X 10 13 Hz < v < few X 10 14 Hz 
marked by Au <C v in the figure, rapid oscillations of the polar- 
ization position angle may not be resolved, thus suppressing linear 
polarization, rendering circular polarization dominant. The calcu- 
lations leading to this result correspond to a typical GRB jet with 
an opening angle 8j 3> T -1 , implying / 2> h (sec Sl3*t. 



over / (see §03) ■ We repeated the calculations of plasma 
propagation effects on the polarization properties of early 
afterglow emission, assuming uniform, close to equipar- 
tition fields (es = 0.1) permeate the reverse and forward 
shocked plasmas, but this time neglecting / in compar- 
ison to h [i.e., setting $ = tt/2 in equations - ©]■ 
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Fig. 6. — Same as Fig. |2J but with i9 set to n/2, implying 
effectively h 3> /■ This situation is a-typical, and is probable for a 



narrow jet of an opening angle 9j 



with the jet observed off- 



axis, and the uniform magnetic field directed normal to the shock 
plane. Plasma parameters correspond to a fireball expanding into 
ISM. However, we find that the picture is qualitatively similar also 
in the wind scenario. 



The resulting linear, circular and total degrees of polar- 
ization of the combined forward-reverse shock emission 
are shown in Fig. |S| The emerging picture is markedly 
different than the one described above (§0J Fig. 0. Lin- 
ear polarization is not suppressed by Faraday rotation, 
and dominates over circular polarization in the entire fre- 
quency range of interest. (Notice, however, the compa- 
rable degree of circular polarization in the emission from 
the forward shock at low frequencies, v < fewx 100 MHz, 
due entirely to the intrinsic circular polarization of syn- 
chrotron at low frequencies) . Additionally, we notice the 
absence of Faraday oscillations of the polarization posi- 
tion angle x a t high (IR to UV) frequencies. Not being 
ironed out by the oscillations, linear polarization is there- 
fore expected to prevail in this case in the IR - UV range 
as well, in contrast to the situation in the case of a typ- 
ical jet. Interestingly, linear polarization is suppressed 
in the microwave and IR band (contributed solely by re- 
verse shock emission) not by the action of plasma effects, 
but rather by the suppression of the Stokes parameter Q 
relative to I at large optical depths 6 . 

4.4. Implications of magnetic field strength and 
structure 

As mentioned above, the fraction of the post-shock en- 
ergy density carried by the magnetic field, es, is one 
of the model parameters least constrained by observa- 
tions. Since the propagation effects addressed in this 
paper are particularly sensitive to the strength of the 
magnetic field, it is natural to explore how the polariza- 
tion of emitted radiation depends on the amplitude of B. 
We therefore consider es = 10~ 4 as a fiducial value for 
the weak field regime, keeping all other fireball parame- 
ters unchanged. 

We find that the general features exhibited by the lin- 
ear, circular and total degrees of polarization are broadly 
similar to those described above for the uniform equipar- 
tition field scenario described above. In particular, we 

6 See, e.g., Pacholczyk (1977) for an analogous suppression of 
U L at v > v m 
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find that the linear polarization is suppressed at low fre- 
quencies. There are some marked differences, though : 
the transition frequency from circular to linear polariza- 
tion is decreased by approximately an order of magnitude 
compared to the strong field scenarios described above. 
Also, the oscillation pattern of the position angle in the 
reverse shock scenario does not survive the decrease in 
B, and is absent in the weak field regime. 

In the case where the structure of the magnetic field 
is entangled over length scales £ much shorter than the 
typical width of the emitting slabs W, the fluctuating 
phase accumulated between the two normal modes av- 
erages to zero, and its r.m.s. is suppressed by a factor 
(i/W) 1 / 2 -C 1, resulting in the absence of the propaga- 
tion effects discussed above. It is helpful in this context 
to separate the discussion of the two fireball and obser- 
vation geometric regimes. Addressing first the typical jet 
regime (0j ^ T" 1 ) in the random field scenario {£ <§; W), 
we conclude that linear polarization is not subject to 
suppression due to the Faraday rotation. Nevertheless, 
this does not imply a high degree of linear polarization, 
because the overall orientation of synchrotron radiation 
emitted from the many different patches averages to zero. 
If, however, the scale of the coherent patch is limited only 
by causality and is thus comparable to the width of the 
emitting slab, viz. £ ~ W (Gruzinov & Waxman 1999), 
Faraday rotation is not completely suppressed, and some 
suppression of 11^ with respect to He is expected. The 
frequencies of transition from circular to linear polariza- 
tion are not expected to change substantially, although 
the degrees of circular and linear polarization will be 
lower by a factor of a few (with a maximal II at the 
level of ~ 10%), due to the different directions that the 
magnetic field may assume in the ~ 50 coherent, causal 
patches apparent to the observer (Gruzinov & Waxman 
1999). Turning our attention to the narrow jet regime 
now, a random, £ <ti W field will suppress all propagation 
effects. In this case a high degree of linear polarization 
— Ul ~ few x 10% — may arise, if the random magnetic 
field has a dominant (fluctuating) component either nor- 
mal to the plane of the shock or parallel to it (Waxman 
2003b; Nakar et al. 2003). The high polarization stems 
from the orthogonality of the direction k of observed ra- 
diation to the direction of the magnetic field in the frame 
comoving with the plasma (see Fig. right for illustra- 
tion). Independent patches limited by causality, like in 
Gruzinov & Waxman (1999), will not suppress propaga- 
tion effects altogether, and the polarization pattern of 
the emitted radiation is expected to follow qualitatively 
either Fig. [SJor Fig. [ft] depending on the dominant com- 
ponent of the fluctuating magnetic field, with maximum 
levels of Hl and He decreased by a factor of few. These 
predictions are summarized in Tabled 

To conclude, then, the polarization of the early after- 
glow emission provides a robust measure of the field's 
complexity : A dominant, high degree of circular polar- 
ization at v < 1 GHz and again at few x 10 10 Hz < v < 
few x 10 14 Hz is a clear signature of a large scale, coher- 
ent field present in the plasma. On the other hand, a 
high degree of linear polarization may result with high 
probability in the narrow jet regime both from a large 
scale field normal to the shock plane, or from an entan- 
gled field, and is therefore less informative as far as field 



structure is concerned. 

5. CONCLUSION 

Ever since it was realized that the dominant mecha- 
nism responsible for the observed radiation from GRB 
afterglows is synchrotron emission, the spectra and po- 
larization at various stages of fireball expansion were cal- 
culated by integration of the relevant emissivities and ab- 
sorption coefficients over varying conditions in the fire- 
ball plasma. Indeed, lightcurves and spectra are suc- 
cessfully reproduced by such procedures. However, the 
propagation effects of a magnetized, relativistic plasma 
on the polarized synchrotron radiation may have sig- 
nificant observable outcomes, and were not treated so 
far. In particular, the bircfringent nature of the magne- 
tized plasma, leading to a phase shift between the normal 
propagation modes, which is accumulated on the coher- 
ence length scale of the magnetic field, and manifested as 
the Faraday effect, may bear significant consequences for 
the polarization of GRB afterglows; A full treatment of 
this problem, taking into consideration the central part 
played by synchrotron cooling (see below), is presented 
here for the first time. 

In this work we generalized results by Sazonov (1969), 
and calculated the terms of the dielectric susceptibility 
tensor responsible for the Faraday effect for a plasma 
with an isotropic, relativistic distribution of electrons. 
We also derived new simple, analytic expressions for the 
propagation coefficients in various frequency ranges (see 
Eqs. 0, [E], and Appendix I A. 1(1 . Similar expressions 
were derived independently by Matsumiya & Ioka (2003). 

We have applied these results to plasma conditions 
typical of the early afterglow, where plasma effects on 
the polarization of propagating synchrotron radiation are 
significant. We considered both the forward and reverse 
shocks, and showed that if the plasma is permeated by a 
large-scale, close to equipartition magnetic field, the po- 
larization pattern of emitted radiation differs decisively 
from the naive estimates usually quoted in the literature. 
We showed that for a typical GRB jet, in which 9j ^> 
during the early afterglow, a large-scale field results in a 
strong suppression of linear polarization at low frequen- 
cies, i.e. v < 3 GHz (y < 3 x 10 13 Hz) in a forward (re- 
verse) shock due to Faraday rotation, and a prevalence of 
circular polarization. Above these transition frequencies 
linear polarization takes over, and gradually approaches 
the canonical value. At the frequencies of transition, the 
total degree of polarization obtains a global minimum 
of 10-20%. The emission from the reverse shock above 
the characteristic synchrotron frequency exhibits strong, 
fast oscillations of the polarization position angle as a 
function of frequency. The oscillations mark the emis- 
sion up to UV frequencies (10 15 — 10 16 Hz); since they 
may be spectrally resolved only above ~ few x 10 14 Hz, 
this implies a suppression of linear polarization (and a 
dominating circular polarization as a result) up to the 
UV band. These prominent fingerprints of the presence 
of a large-scale magnetic field in the plasma are rather 
independent of the environment into which the fireball 
expands, whether it has a uniform density (ISM) or a 
wind density profile (see § 14.2(1 . The suppression of linear 
polarization at low frequencies was shown to characterize 
also the case of a uniform field which is much below the 
equipartition value. However, in the weak field regime 
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the frequency of transition from circular to linear polar- 
ization decreases by approximately an order of magni- 
tude, and the oscillation pattern of the position angle is 
absent (whence also the dominance of 11c up to the UV 
band). These propagation effects are expected to vanish 
if the field is entangled over length scales much smaller 
than the extent of the emitting plasma. 

A large-scale magnetic field is likely to give rise to dif- 
ferent dispersive effects on the polarization properties if 
the jet is narrow (6j ~ T^ 1 ) and observed off-axis, and 
if additionally the field is predominantly normal to the 
shock plane. In this case we predict dominance of linear 
polarization at all frequencies (with the exception of a 
comparable circular component at low radio frequencies). 
However, this is not an exclusive hallmark of a large scale 
magnetic field in the narrow, off-axis jet regime, since a 
high degree of linear polarization may result also if the 
field is entangled over small length scales, as may have 
been the case in GRB 021206. 

Recently, Matsumiya and Ioka (2003) treated similar 
propagation effects in the context of GRB afterglows, 
neglecting the inhomogeneity of the emitting slabs due 
to synchrotron losses. As we demonstrated explicitly in 
§ 14.21 cooling has a decisive effect, most importantly on 
the values of the propagation coefficients. This leads 
to differences between Matsumiya and Ioka's results and 



ours (e.g. the higher degree of circular polarization, by 
approximately an order of magnitude, that we find in the 
transition frequencies 1 GHz and few xlO 13 Hz). 

Our results suggest that spectro-polarimetric observa- 
tions in the radio and IR bands during the early stages of 
the afterglow may be used as a unique probe on the struc- 
ture and strength of the magnetic field. A handle on the 
field configuration is of particular importance in study- 
ing the plasma conditions in the reverse shock, since this 
plasma is essentially the one outflowing from the com- 
pact source, and hence may reveal valuable information 
on the processes driving the outflow and on the nature of 
the progenitor. This technique, complementary to mea- 
surements of the 7— ray polarization, may therefore be 
powerful in placing stringent constraints on models of 
the "inner engine" and the origin of the magnetic field, 
and promote our understanding of the physics of colli- 
sionless shock waves. 
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APPENDIX 

PROPAGATION, EMISSION AND ABSORPTION COEFFICIENTS IN A RELATIVISTIC PLASMA 

The propagation coefficients in a relativistic plasma 

We briefly review here the derivation of the coefficients / and h [see Eqs. - ||I2J], associated with propagation 
effects, for the case of a plasma with an isotropic, relativistic electron population. The deviation of the dielectric tensor 



f(",k) 



47TKy(u;,k) 



(Al) 



from its value in vacuum (8ij) is assumed to be small, so that the radiation can be regarded as transverse (Sazonov 
1969). 

To obtain /, a quantity which is first order in B, we separate the perturbation to the equilibrium distribution function 
into two contributions, one which is independent of the magnetic field, and one which is linear in the magnetic field. 
The Maxwell equations lead to a dispersion relation : 

c 2 k 2 /cu 2 E l =e l3 E 3 , z,j = l,2, (A2) 
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which, when coupled to a linearized Vlasov equation, gives the non-diagonal elements : 



ei2 



-621 



cosz? 



dV 



1 



d/(°V ±1 



(A3) 



uj J (lo — fcu||) 2 dp "f 2 p 

,, — Amn e e 2 jm e and ujb = eB/m e c are the non-relativistic (electron) plasma frequency and Larmor frequency, 



Here uj 2 

respectively; ■& is the angle between k and B, 7 = y/p 2 + m 2 c? /m e c is the particle's Lorentz factor, and f^°'(p) is 
the unperturbed isotropic electron distribution function. 

Sagiv and Waxman (2002) showed numerically that the value of 612 is not very sensitive to the exact shape of 
the distribution function, and may be approximated by substituting a mono-energetic electron distribution /*- -'(7) oc 
(5(7 — 7 m ) in equation (|A3|I . with j m a characteristic electron Lorentz factor (typically at the peak of the distribution 
function). This leads to a simple expression for the propagation coefficient / in the rclativistic regime : 

e 3 n e .B cos 1? In 7 m 

Uei(v) - yITIl — ' A4 ) 

To obtain h re i, we started with equation (9) of Sazonov (1969), who used a kinetic equation approach to obtain 
analytic expressions for the elements of the susceptibility tensor «y. Using the notation of the current paper, this 
equation reads : 



1 e 2 m 2 c® 

1 U) 2 



lub sin i9 



d£ 



d£ 



dn(£)/d£ 



£ 2 



l A Y{z) , (A5) 
2//3 is a normalized 



with £ = ^m e c the particle energy, and dn(£)/d£ the electron energy spectrum; z = {jv/v syn 

frequency, where v syn = 3eBj 2 / 4irm e c is the characteristic synchrotron frequency for electrons with Lorentz factor 7; 

The function Y(z) is defined as : 

z 2 f C°° dv 
Y { z ) = ^{ I —exp(ivz + w 3 /3) 



3^ 
2 



/OO 
dv v exp(iisz + zi/ 3 /3) 
-00 



dv v cos(i^z -(- ^ 3 /3) , 



(A6) 



where V denotes Cauchy's principal value of the integral. 

Since we are interested in h, only the real part of Y{z) matters, and it can be easily shown that in fact only the last 
integral contributes. At frequencies v <C v mi Ku is dominated by electrons at the peak of the distribution function. 
We therefore take dn(£)/d£ oc 5(£ — j m m e c 2 ). Expanding i?e[Y"(z)] at the limit z <C 1, we finally obtain 

eVS 2 sin 2 i? 7 „ / << N 4/3 



[y < v m ) 



(A7) 



2 \ Vm (^77 

At frequencies v > v m the value of h re i is sensitive to the shape of the electron distribution function. Specifically, 
for a power-law distribution, we simplify equation (20) of Sazonov (1969) to get the second of equations I© : 



e A n P B 2 sin 2 d 



In 



(p-2)/2 



(v > V m ) 



(A8) 



47T 2 TO 3 C 3 ^ 3 \p — 2 J "" \V„ 

The emissivity and absorption coefficients for synchrotron radiation 

For completeness we provide expressions for the emissivity and absorption coefficient of the synchrotron radiation, 
using the notation x — v/v syn [see, e.g. Pacholczyk (1977), but note the typos in this source]. 

( £I \ 



V3 e 3 B 



47r m e c 2 



sin i3 



dn(£ e 



VvJ 



2cottf 



V 7e 

VSe 3 B 



d£ e 

x XT K 5/z{z)Az 
xK 2 /a(x) 



sinz9— ^ 



K 1/3 {z)dz 
<1 fdn{£ e )/d£, 



d£. 



(A9) 



87r m e v- j 
' x JT K 5/3(z)dz 



d£ ( 



£' 



d£ e 



(A10) 



Note that for a synchrotron source zjj = 0, u = [Pacholczyk, (1977)], and we may neglect v in comparison to k 
and q. 



